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Abstract 

We present the results from a large near-infrared spectroscopic survey with Subaru/FMOS {FastSound) consist¬ 
ing of ^ 4,000 galaxies at z ~ 1.4 with significant Ha detection. We measure the gas-phase metallicity from the 
[N ll]A6583/Ha emission line ratio of the composite spectra in various stellar mass and star-formation rate bins. 

The resulting mass-metallicity relation generally agrees with previous studies obtained in a similar redshift range 
to that of our sample. No clear dependence of the mass-metallicity relation with star-formation rate is found. Our 
result at z ^ 1.4 is roughly in agreement with the fundamental metallicity relation at z ~ 0.1 with fiber aperture 
corrected star-formation rate. We detect significant [S ll]AA6716,6731 emission lines from the composite spectra. 

The electron density estimated from the [S ll]AA6716,6731 line ratio ranges from 10 - 500 cm“^, which gener¬ 
ally agrees with that of local galaxies. On the other hand, the distribution of our sample on [N ll]A6583/Ha vs. 

[S ll]AA6716,6731/Ha is different from that found locally. We estimate the nitrogen-to-oxygen abundance ratio 
(N/O) from the N2S2 index, and find that the N/O in galaxies at z ^ 1.4 is significantly higher than the local values 
at a fixed metallicity and stellar mass. The metallicity at z 1.4 recalculated with this N/O enhancement taken into 
account decreases by 0.1 - 0.2 dex. The resulting metallicity is lower than the local fundamental metallicity relation. 

Keywords: galaxies: evolution — galaxies: abundances — galaxies: high-redshift — galaxies: ISM 


1. Introduction 

The metallicity in the gas-phase (hereafter, metallicity) is 
one of the key parameters required for understanding galaxy 
formation and evolution. Heavy elements produced in stars 
through the star-formation activity are returned into the inter¬ 
stellar medium (ISM), which is also regulated by gas infall and 
outflows. Thus, the metallicity can be a tracer of the past star- 
formation history including the gas flow history (e.g.,Tinsley 
1980). Enormous efforts on studying the correlation between 
the stellar mass and metallicity (hereafter mass-metallicity re¬ 
lation) over a wide dynamic range has been made at z ^ 0.3 
(e.g., Fequeux et al. 1979, Tremonti et al. 2004, Gallazzi 2005, 
Fee et al. 2006, Kewley & Ellison 2008, and Andrews &. 
Martini 2013), at intermediate redshift of 0.3 ^ z ^ 1 (e.g., 

* Based on data collected at Subaru Telescope, which is operated by the 

National Astronomical Observatoiy of Japan. 


Savaglio et al. 2005, Rodrigues et al. 2008, Famareille et al. 
2009, Zahid et al. 2013, Fy et al. 2014), and at higher red- 
shifts, z ^ 1 (e.g., Erb et al. 2006, Maiolino 2008, Mannucci 
et al. 2009, Hayashi et al. 2009, Yabe et al. 2012, Belli et 
al. 2013, Henry et al. 2013, Yabe et al. 2014, Zahid et al. 
2014b, Cullen et al. 2014, Steidel et al. 2014, and Sanders et al. 
2015). The compilation of these studies shows that the mass- 
metallicity relation evolves with redshift such that the metallic¬ 
ity increases gradually with decreasing redshift at a fixed stellar 
mass (Maiolino 2008, and Zahid et al. 2013). 

The intrinsic scatter of the mass-metallicity relation and its 
dependence of the star-formation activity has also been inves¬ 
tigated (e.g., Tremonti et al. 2004, Ellison et al. 2008, Lara- 
Lopez et al. 2010b, Mannucci et al. 2010, Yates, Kauffmann, 
& Guo 2012, and Andrews & Martini 2013). A fundamental re¬ 
lation between stellar mass, metallicity, and star-formation rate 
(SFR) for Sloan Digital Sky Survey (SDSS) galaxies at z 0.1 
has been proposed by Mannucci et al. (2010), with an exten- 
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sion to galaxies of lower stellar mass (Mannucci, Salvaterra, 
& Campisi 2011, Hunt et al. 2012, and Zahid et al. 2012). 
According to this fundamental metallicity relation (hereafter, 
FMR), galaxies with higher SFR show lower metallicity at a 
fixed stellar mass. The evolution of the mass-metallicity rela¬ 
tion can be explained by the evolution of the SFR of galaxies 
(Mannucci et al. 2010). The presence of a FMR and its origin 
has, however, been argued in many studies (e.g., Mannucci et 
al. 2010, Bothwell et al. 2013, Forbes et al. 2014, and Zahid et 
al. 2014a). 

Recently, the presence of the FMR at higher redshift has 
been explored in various studies (e.g., Cresci et al. 2012, Yabe 
et al. 2012, Vergani et al. 2012, Wuyts et al. 2012, Stott et al. 

2013, Yabe et al. 2014, Zahid et al. 2014b, Stott et al. 2014, 
Divoy et al. 2014, and Steidel et al. 2014). The dependence of 
the mass-metallicity relation on the SFR at high redshift, how¬ 
ever, still remains unclear. Some studies suggest a dependence 
of the mass-metallicity relation on SFR or specific SFR within 
their samples (Yabe et al. 2012, Stott et al. 2013, Stott et al. 

2014, and Zahid et al. 2014b), but others do not show such a 
clear dependence (Yabe et al. 2014, Steidel et al. 2014, and 
Sanders et al. 2015). This disagreement may be partly due to 
the smallness of the samples, as the number of galaxies in each 
study is still limited up to ^300. Furthermore, no clear SFR 
dependence is due to the narrow range of the SFR (see discus¬ 
sions by Stott et al. 2013; Yabe et al. 2014). Furthermore, some 
studies suggest offsets at high redshift from the local FMR. By 
using their sample at z = 0.84 - 1.6, Yabe et al. (2012), Stott 
et al. (2013), and Yabe et al. (2014) reported that the metallic¬ 
ity is generally higher than the extrapolation of the local FMR 
to the lower values of p (= log(M*) — alog(SFIi)), where a 
is the projection parameter of the FMR (see Mannucci et al. 
2010). On the other hand, Cullen et al. (2014) also reported the 
deviation from the local FMR at z ^ 2 — 4 with R23 metallic¬ 
ity indicator, but they suggest a counter trend; the metallicity 
is lower than the local relation. A Larger sample spanning a 
wider range of parameter space is needed to reveal the FMR 
and other parameter dependence of the mass-metallicity rela¬ 
tion at high redshift. 

The FastSoimd project (Tonegawa et al. 2015b) is a cosmo¬ 
logical redshift survey to measure redshifts of several thou¬ 
sand star-forming galaxies at z = 1.18 - 1.54 in the total area 
of ~20 deg^ in the Canada-France-Hawaii Telescope Legacy 
Survey (CFHTLS) wide fields by using the Ha emission lines 
obtained with Fiber Multi-Object Spectrograph (FMOS) on the 
Subaru Telescope. Although the main science goal of the sur¬ 
vey is to measure the redshift space distortion in this redshift 
range, the survey also provides for ancillary science cases with 
a large near-infrared (NIR) spectroscopic sample of Ha emit¬ 
ting galaxies at z > 1. The targets in this survey are galaxies 
which are relatively bright in Ha emission line, observed in 
a short exposure time, i.e., the target galaxies have moderately 
high SFR. In this work, by using an unprecedentedly large NIR 
spectroscopic sample, we examine the mass-metallicity rela¬ 
tion of star-forming galaxies and explore the high-SFR end of 
the FMR. 

Throughout this paper, we use the following cosmological 
parameters: Hq = 70km s“^ Mpc“^, Ha = 0.7, and = 0.3. 
All magnitudes given in this paper are in the AB magnitude 
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Fig. 1. The distribution of 2 :'-band magnitude and g' — r' color for 
target sample (blue contour), Ha detected sample (red contour), and 
the K-selected photometric redshift sample at 1.2 < Zph < 1.5 in the 
SXDSAJDS field (black contour). The contour levels are 68%, 90%, 
and 97%, respectively. In the top and right panels, the peak-normal¬ 
ized histograms of the magnitude and color are shown for each sample. 

system. We use the Salpeter initial mass function (IMF). 

2. Sample Selection and Observations 

The target sample in the FastSound survey is selected by 
photometric redshift and the predicted Ha flux based on the 
SFR and color excess from SED fitting with u*, g', r', i', and 
z' photometric data of the four CFHTLS wide fields (Wl, W2, 
W3, and W4). The photometric redshifts of the galaxies are 
derived using LePhare (Ilbert et al. 2006). As a target sam¬ 
ple, galaxies with 1.1 < Zpt < 1-6, 20.0 < z' < 23.0 mag, and 
g' — r' < 0.55 mag are selected. Furthermore, an additional cri¬ 
terion of an expected Ha flux of ^ 1.0 x 10“^® erg s“^ cm“^ 
is applied. The magnitude, color, and the expected Ha flux cut 
for the target selection are made aiming to enhance the detec¬ 
tion rate in the survey. It should be noted that some selection 
biases are raised by the expected Ha flux cut, though the effect 
is marginal according to the previous studies (Yabe et al. 2012, 
and Yabe et al. 2014). The detailed descriptions on the sample 
selection are presented by Tonegawa et al. (2015b). In figure 1, 
the distribution of our target sample in z'-band magnitude vs. 
g' — r' color diagram is shown with more completed sample at 
1.2 < Zph < 1.5 in the SXDS/UDS held by Yabe et al. (2014). 
The peak-normalized histograms show that our target sample 
is relatively biased toward brighter magnitude and bluer color 
compared to the SXDS/UDS sample, because of our magni¬ 
tude and color selection. Dusty and/or old galaxies may be left 
out of our sample selection. 

The observations were carried out in 40 nights from 
September 2011 to January 2014 by using Subaru/FMOS, 
which is a fiber-fed multi-object NIR spectrograph with 400 
fibers (1.2" diameter) at the prime focus of the Subaru 
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Telescope (30' diameter). Details of the instrument are pre¬ 
sented by Kimura et al. (2010). We use the spectrograph set¬ 
ting of the high resolution mode (H-short prime; R ^ 2400), 
covering the wavelength range from 1.45 pm to 1.67 pm. We 
used the normal beam switching mode during the observations, 
in which all the fibers are available for the target objects (up 
to ^ 400). In the typical observational process of this survey, 
after the initial fiber configuration, a set of on-target and off- 
target frame with exposure time of 15 minutes each is taken. 
The typical beam-switch offset is 10" - 15" from the target. 
After re-configuration of the fiber positions, another set of the 
on- and off-frame is taken in the same manner. The total on- 
source exposure time is 30 minutes for each object. The seeing 
size, which was monitored with the optical sky camera during 
the observations, ranged from ^ 0.5" to ^ 1.7" with a median 
value of 1.08". We observed 118 field-of-views in the four 
CFHTLS wide fields (^20 deg^), covering 38,615 targets in 
total. 

The basic data reductions are carried out by using the stan¬ 
dard FMOS pipeline FIBRE-pac (Iwamuro et al. 2012) '. The 
emission line detection from the extracted ID spectra is carried 
out by using an automated emission line detection software 
(FIELD), which was developed and tested for the FastSound 
survey (Tonegawa et al. 2015a). We detect Ha emission 
lines from 9,803 galaxies with signal-to-noise ratio (S/N)>3.0, 
which include 4,139 and 2,740 galaxies with S/N>4.0 and 5.0, 
respectively. According to the analysis by Tonegawa et al. 
(2014) using inverted images, the contamination rate due to 
false detections is ~ 5% for a S/N cut of 4.0. In this work, we 
use the S/N threshold of 4.0 for Ha detection. The spectro¬ 
scopic redshifts of the detected objects range from 1.18 to 1.54 
with the average (median) value of 1.362 (1.356). The redshift 
distribution is presented by Okada et al. (2015). The observed 
line width (FWHM”*’'*) ranges from ^ 100 to ^ 1000 km s“^ 
with the average (median) value of 232 (212) km s“^. After 
deconvolving the spectral resolution (R ^ 2400) by subtract¬ 
ing in quadrature the FWHM measured from unresolved sky 
lines, the average (median) intrinsic line width (FWHM'"'^) is 
203 (187) km s“^. Figure 1 shows that the Ha detected sample 
has similar z'-band magnitude and g' — r' color range to those 
of the target sample. 

Since our observations targeting the Ha detection do not 
cover the J-band, we have no information on H/3 and 
[O III] AA4959,5007. The emission line ratio diagnostics us¬ 
ing [N II]A6583/Ha and [O lll]A5007/H/3 (Baldwin, Phillips, 
& Terlevich 1981) and the recent mass-excitation diagnostics 
using stellar mass and [O III] A5007/H/3 line ratio (Juneau et al. 
2014) cannot be used. According to diagnostics by Kewley 
et al. (2006), in the range of the obtained [N ll]A6583/Ha 
in this work, a spectrum is identified as a star-forming re¬ 
gion if log([0 lll]A5007/H/3) is ^ 0.7 — 0.8. In the stacked 
spectrum of the whole sample, which will be described in de¬ 
tail in the next section, [O l]A6300 emission line is detected 
with S/N^20. According to Kewley et al. (2006), a measured 
log([0 l]A6300/Ha) of -1.55 is identified with star-forming 

* Since the wavelength calibration of the FMOS data is carried out with vac- 

uum wavelength, we use the vacuum wavelength of each emission line in 

the analysis of this work. Thi'oughout this paper, however, we note the 

wavelength of each line in air, such as [N II]A6583, [S II]AA6716,6731, etc. 
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Stellar Mass [Mq] (without NIR/MIR) 


Fig. 2. Comparison between stellar mass derived with NIR/MIR pho¬ 
tometry and that without NIR/MIR photometry. Red filled cir¬ 
cles are objects with UKIDSS/DXS data and blue filled triangles 
are objects with UKIDSS/UDS data. Dotted line shows an equiv¬ 
alent line and dashed line indicates the linear regression line. 

galaxies if log([0 lll]A5007/H^) is ^ 0.6. In the study by 
Yabe et al. (2014), most of the sample galaxies in the red- 
shift range of our sample show log([0 lll]A5007/H/3);< 0.6. 
In any event, active galactic nuclei (AGNs) with these low 
[N ll]A6583/Ha and [O l]A6300/Ha are very rare in the local 
universe. Furthermore, as we will discuss in subsection 5.1, the 
[N II] A6583/Ha and [S ll]AA6716,6731/Ha line ratios indicate 
that our sample does not suffer severely from AGN contamina¬ 
tion. As we mentioned above, most of the line widths of our 
sample are ^1000 km s“^, which means that our sample is not 
severely contaminated by broad line AGNs. In this work, we 
exclude only a marginal fraction (^1.1 %) of broad line emis¬ 
sion galaxies (FWHM°^® > 1000 km s“^) from the analysis 
below. Thus, the total number of galaxies used in this work 
after excluding the possible broad AGN candidates is 4,094. 

3. Measurements of stellar mass, star-formation rate, and 
metallicity 

3.1. Stellar masses 

The stellar masses used in the present analysis are derived 
from multi-wavelength spectral energy distribution (SED) fit¬ 
ting using photometric data from the rest-frame ultraviolet 
(UV) to mid-infrared (MIR). The initial catalogue of our 
sample with optical photometry is first cross-matched to the 
UKIRT (United Kingdom Infrared Telescope) Infrared Deep 
Sky Survey (UKIDSS) Deep Extragalactic Survey (DXS) cat¬ 
alogue derived from UKIRT Wide Eield Camera (WECAM) 
imaging data in the J- and AT-band near-infrared (NIR) bands 
(Lawrence et al. 2007), which covers part of CEHTLS W1 and 
W4. In total, 487 objects out of 701 objects are detected in the 
J and/or K bands. The total magnitude of the detected objects 
in the NIR with S/N>5 ranges from J ~ 20.0 - 23.0 mag and 
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Fig. 3. The distribution of the stellar mass and the SFR of our 
sample. Contours for 96%, 85%, 66%, 24%, and 4% level are 
shown in gray cuiwes. The average value in each stellai' mass 
and SFR bin (indicated by horizontal and vertical dashed lines) 
is shown by a red circle. The number of galaxies in each bin 
is also presented at the bottom left comer in each section. The 
dotted line is the galaxy main sequence at z ^ 1.4 by Whitaker 
et al. (2012). In the top and right panels, the peak-normalized 
distribution of stellar mass and SFR is shown, respectively. 

K ~ 19.5 - 22.8 mag. Our sample is also cross-matched to 
the Spitzer Enhanced Imaging Products^ (SHIP) catalogue de¬ 
rived from Spitzer Infrared Array Camera (IRAC) 3.6pm (ch. I) 
and 4.5pm (ch.2) MIR data, where combined images are con¬ 
structed from several earlier programs. In total, 422 objects of 
our sample are detected in ch.l and/or ch.2 band. The total 
magnitude of the detected objects in MIR with S/N>5 ranges 
from ch.l ~ 19.0 - 22.6 mag and ch.2 ^ 18.8 - 22.8 mag. 791 
objects (~ 20% of the total sample) are detected in at least one 
of the NIR to MIR bands. 

SED fitting was undertaken using the SEDfit code (Sawicki 
2012) as described by Yabe et al. (2014). Here, we assume 
the Salpeter IME (0.1 - 100 Mq) and the Calzetti extinction 
law. Eor objects without NIR or MIR data, stellar masses are 
estimated using only optical data from the u*- to z'-bands. 
To evaluate the uncertainties arising from the lack of NIR 
and/or MIR (hereafter, simply NIR/MIR) data, we compare 
stellar masses derived with and without NIR/MIR data by us¬ 
ing the subsample with both optical and NIR/MIR photome¬ 
try in figure 2. Although the stellar masses derived without 
using the NIR/MIR data are larger by up to ^ 1 dex in the 
most massive end, the agreement is reasonably good in the 
mass range from ^ 10® to a few x 10^® Mq, in which most 
of our sample galaxies are. We discuss the possible effect 
of the uncertainty in the stellar mass derivation on results in 
section 4.1. Using a linear function to correct this bias, we 
find; = 0.63 x -f 3.50, 

dashed line in figure 2. with a scatter corresponding to cr = 0.22 

^ http;//irsa.ipac.caltech.edu/data/SPITZER/Enhanced/SEIP/overview.html 


dex. In the following analysis, if the stellar mass can only be 
derived from optical data, we correct it applying this relation. 

The comparison of the stellar mass derived with and with¬ 
out NIR/MIR data in figure 2 could be biased if the nature of 
the SED of the NIR/MIR undetected sample differs from that 
where NIR/MIR detection is available. This bias could be crit¬ 
ical for the less massive portion of figure 2. To address this, 
we check the stellar masses for a subset of our sample in the 
UKIDSS Ultra Deep Survey (UDS DR8; Lawrence et al. 2007) 
field, where deeper NIR data is available. The UDS field cov¬ 
ers part of CEHTLS W1 but is not overlapped with DXS field. 
In this deeper NIR sample, 115 out of 140 of our objects are 
detected down to ^23.6 mag with S/N^IO. The stellar masses 
derived with the deeper NIR data are compared to that without 
the NIR data. Eigure 2 shows that the comparison using deeper 
UDS data remains consistent with the regression line derived 
using the DXS data. 

3.2. Star-formation rates 

The star formation rate (SER) is estimated from the extinc¬ 
tion corrected Ha luminosity according to the methodology 
discussed by Yabe et al. (2012) and Yabe et al. (2014). The 
color excess for the stellar continuum is derived from u* — r' 
color, corresponding to the rest-frame UV for the redshift range 
of our sample, with an extinction law by Calzetti et al. (2000). 
Based on the prescription by Cid Eernandes et al. (2005), we 
correct the broad-band color excess to that for nebular emis¬ 
sion, applying a factor of ^ 2. We use the relation by Kennicutt 
(1998) to infer the SER from the Ha luminosity. Eollowing the 
study by Tonegawa et al. (2014) using a sample similar to the 
present one, we applied a fiber aperture correction of 0.5. The 
uncertainties of our flux estimates will be described in Okada et 
al. (2015). The SERs derived from our Ha luminosities agree 
with those inferred from the dust-corrected rest-frame UV lu¬ 
minosity density with a scatter of only cr = 0.26 dex. The distri¬ 
bution of stellar masses and SERs of our sample is presented in 
figure 3 and summarized in table 1. It is worth noting that our 
sample mostly comprises galaxies with higher SERs compared 
to the ‘main sequence’ at z ^ 1.4 (Whitaker et al. 2012). 

3.3. Metallicities 

The metallicities in our analysis are derived from the 
[N Il]/Ha line ratio according the relation introduced by Pettini 
& Pagel (2004); 

12H-log(0/H) = 8.90-f 0.57 xN2, (1) 

where N2 = log([N ll]A6583/Ha). Although it is known that 
the N2 index saturates above solar metallicity, most of our sam¬ 
ple has sub-solar metallicity and the effect of this saturation 
will be small. In our analysis, the errors do not include any pos¬ 
sible intrinsic error in the calibration (~0.2 dex; Pettini & Pagel 
2004) or potential systematic differences between the vari¬ 
ous metallicity indicators (up to ^0.7 dex; Kewley & Ellison 
2008). 

[N ll]A6583 emission lines from 224 objects are detected 
with S/N > 3, while a large part of our sample shows no de¬ 
tection. In order to measure the average metallicity and estab¬ 
lish the mass-metallicity relation, we apply a spectral stacking 
analysis follow the method described by Yabe et al. (2012) and 
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Table 1. Stellar mass and SFR of each bin in our sample. 


Bin 

Num.“ 

log(M*)'> 

[Mq] 

log(SFR)'’ 
[Mq yr-1] 

1 

171 

9.55 

1.81 

2 

208 

9.90 

1.81 

3 

167 

10.05 

1.82 

4 

118 

10.26 

1.83 

5 

148 

10.61 

1.81 

6 

181 

9.57 

2.05 

7 

187 

9.89 

2.03 

8 

158 

10.06 

2.03 

9 

150 

10.24 

2.04 

10 

134 

10.53 

2.03 

11 

205 

9.61 

2.21 

12 

146 

9.90 

2.21 

13 

171 

10.08 

2.21 

14 

156 

10.23 

2.21 

15 

136 

10.55 

2.21 

16 

154 

9.61 

2.38 

17 

137 

9.89 

2.38 

18 

168 

10.07 

2.37 

19 

190 

10.25 

2.39 

20 

162 

10.52 

2.39 

21 

101 

9.65 

2.62 

22 

133 

9.86 

2.71 

23 

149 

10.08 

2.65 

24 

197 

10.25 

2.65 

25 

235 

10.55 

2.67 


Note - (a) the number of galaxies in each stellar mass and 
SFR bin. (b) median values in each bin. 

Yabe et al. (2014). We construct a weighted-mean composite 
spectrum in each stellar mass and SFR bin. Since the spec¬ 
tra in this study were all obtained in FMOS’ high-resolution 
mode, no correction for the obscuring effect of the OH-masks 
is required. However, some flux may be lost if the line is close 
to a portion of the OH-mask. For Ha and [N ll]A6583, we 
estimated this loss using Gaussian profiles at the appropriate 
emission line wavelength and line width. The average losses 
for Ha and [N ll]A6583 are 7.4% and 13.4%, respectively. The 
effect on the [N ll] A6583/Ha ratio is only 6.5%, corresponding 
to ^ 0.016 dex in 12H-log(0/H), and thus no correction was ap¬ 
plied. The composite spectra in the wavelength range around 
Ha and [N ll]AA6548,6583 lines are presented in figure 4. 

Ha and [N ll]AA6548,6583 line fluxes are measured by spec¬ 
tral fitting of the stacked spectra with multi-Gaussian profile. 
The redshift and line width of other emission lines are based 
on that determined for Ha. We considered various ways of 
spectral stacking (standard mean, median, peak normalized) 
but little difference in the resulting line fluxes. The error on the 
metallicity derived from the stacked spectra is based on boot¬ 
strap resampling with 500 realizations. The metallicity in each 
bin is measured from the obtained stacked spectrum by using 
the N2 method described above, which is also summarized in 
table 2. 


3.4. [S ll]XX6716,6731 emission lines 

The wavelength range of our obtained spectra covers 
[S It] AA6716,6731 lines. The ratio of these two lines is 
sensitive to the electron density (e.g., Osterbrock 1989). 
The composite spectra in the wavelength range around 
[S ll]AA6716,6731 lines are presented in figure 5, respectively. 
We have significant detections in the composite spectra. The 
line fluxes are measured by using the same spectral fitting as 
that for Ha and [N ll] lines, which is also summarized in table 
2. We discuss the electron density in subsection 4.3. 

4. Results 

4.1. The Mass-Metallicity Relation of Star-forming Galaxies 
at z 1.4 

The left panel of figure 6 shows the distribution of the ob¬ 
tained metallicity and stellar mass at ^ ~ 1.4. We also show 
the results from the stacking analysis with different color de¬ 
pending on the SFR. No clear SFR dependence of the mass- 
metallicity relation is seen in the left panel of figure 6. The 
mass-metallicity relation of our sample generally agrees with 
those of previous results with ~340 galaxies at 2 ; ~ 1.4 (Yabe 
et al. 2014) and ~160 galaxies at 2 ; ^ 1.6 (Zahid et al. 2014b), 
where the spectra are taken with FMOS and the metallicity is 
derived with the same N2 method. For this comparison, we 
adopted a Salpeter initial mass function (IMF). In the mas¬ 
sive end with M* ^ 10^°'^ Mq, metallicities of our sample 
are lower than the earlier estimates by up to ^ 0.2 dex. This 
may be partly due to the lack of red galaxies from our sample 
as we mentioned in section 2; galaxies at z ~ 1.4 — 1.6 with 
the redder color and/or the larger color excess show the higher 
metallicity (Yabe et al. 2014; Zahid et al. 2014b). 

In the left panel of figure 6, we compare our result at 2 : ~ 1.4 
to that at 2 : ~ 0.1 obtained with the sample by Tremonti et al. 
(2004) but by using the N2 method (Erb et al. 2006). Here, stel¬ 
lar mass is converted to that obtained assuming Salpeter IMF. 
Our result and previous results by Yabe et al. (2014) and Zahid 
et al. (2014b) at 2 : ~ 1.4 — 1.6 show lower metallicities than the 
result at z ^ 0.1 at a fixed stellar mass by up to ^ 0.2 dex. The 
difference is generally larger in the lower stellar mass, which 
is similar to previous results at z ~ 1.4 — 1.6. 

There exist various uncertainties in the obtained mass- 
metallicity relation at z ^ 1.4. As we mentioned in subsec¬ 
tion 3.1, the stellar mass of a large part of our sample is cor¬ 
rected because of the lack of the NIR/MIR data set. In addition 
to the typical correction factor of ~0.2 dex, there is a scat¬ 
ter of (Tm, = 0.22 dex. There also exists the uncertainty in 
the SFR estimation of <jsfr = 0.26 dex as described in sub¬ 
section 3.2. We estimate the possible effect on the obtained 
mass-metallicity relation arose from these uncertainties. We 
resample randomly on the stellar mass vs. SFR plane in figure 
3 assuming a Gaussian distribution with ctm. and <jsfr men¬ 
tioned above, and then apply the same stacking analysis de¬ 
scribe in subsection 3.3. We repeat this procedure 500 times. 
The standard deviation of the realized metallicities is ^ 0.03 
dex, which is almost similar to the error based on the bootstrap 
method as made in subsection 3.3. 
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Fig. 4. A composite spectrum of our sample at 2 : ~ 1.4 in each stellar mass and SFR bin. The stellar mass increases from left to right 
and the SFR increases from bottom to top. The observed spectra and the best-fit model spectra are shown by red solid lines and blue 
dashed lines, respectively. The positions of Ho and [N II]AA6548,6583 are indicated by vertical dotted lines. The median values of 
log(M*[M q]), log(SFR[MQ yr“^]), and 12-i-log(0/H) are shown at the top-left, top-right, and middle-left part of each panel, respectively. 
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Fig. 5. Same as figure 4, but in the wavelength region around [S II]AA6716,6731 lines, which are indicated 

by vertical dotted lines. The median values of log(M+[M 0 ]), log(SFR[MQ yr“^]), and [S II]AA6716,6731 

line ratio are shown at the top-left, top-right, and middle-right part of each panel, respectively. 
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Fig. 6 . (Left) The metallicity against the stellar mass. Individual object with both Ha and [N II] detection is indicated by a gray dot. The typical errors 
of the stellar mass and metallicity are 0.32 dex and 0.24 dex, respectively. The results from the stacked spectra in 10 stellar mass bins are indicated by 
open circle. The results from the stacked spectra in different SFR (in Mq yr“ ^) bins are shown by filled circles color-coded by the SFR. Filled triangles 
and open squares are results in the past FMOS observations by Yabe et al. (2014) and Zahid et al. (2014b), respectively. For comparison, the result at 
z ~ 0.1 obtained with the sample by Tremonti et al. (2004) by using the N2 method (Erb et al. 2006) is shown as a dashed curve. (Right) The distribution 
of the metallicity against the projected axis fi of the FMR. Here, we use the projection pai'ameter of a = 0.32 for our sample and that by Yabe et al. 
(2014) and Zahid et al. (2014b), which is proposed for SDSS galaxies at z ~ 0.1 by Mannucci et al. (2010). For compailson, the results by Yabe et al. 
(2014) and Zahid et al. (2014b) are also presented as filled triangles and open squares, respectively. The FMR at z ~ 0.1 by Mannucci et al. (2010) is 
presented by thick solid curve. For /r ^ 9.4, the solid curve is extrapolated linearly, which is shown by the thick dashed line. Two FMRs at z ~ 0.1 by 
Zahid et al. (2014b) and Andrews & Martini (2013) are presented by dot-dash curve and dashed line, respectively, where the SFR corrected for the SDSS 
fiber aperture effect is used. We use a projection parameter of a = 0.32 for these two samples. The metallicities are derived by using the N2 method. 


4.2. The Fundamental Metallicity Relation 

The right panel of figure 6 shows the metallicity along the 
projected axis (p = log(M*) — Q!log(SFR)) of the FMR defined 
by Mannucci et al. (2010). Here, we use a projection parame¬ 
ter of a = 0.32, which minimizes the scatter of the local mass- 
metallicity relation. The metallicity of the sample at 2 ; ~ 0.1 
based on the calibration by Maiolino (2008) is converted to 
[N h]A 6583 to Ha line ratio according to the Maiolino (2008) 
calibration, and then is converted to the metallicity again with 
the N2 calibration by Pettini & Pagel (2004). The FMR at 
2 ; ~ 0.1 defined by Mannucci et al. (2010) is shown as thick 
solid curve, and the extension toward lower p, based on a lin¬ 
ear extrapolation, is shown as a dashed line. 

While the SFR used in the FMR by Mannucci et al. (2010) 
is not corrected for the SDSS fiber aperture effect, the SFR 
in our sample is corrected for the FMOS fiber aperture effect 
as described in subsection 3.2. Recently, Andrews & Martini 
(2013) and Zahid et al. (2014b) derived FMRs based on the 
SDSS sample by using SFRs corrected for the aperture effect. 
The resulting FMRs are shown in the right panel of figure 6 
by a dashed line and a dot-dash curve, respectively. Both re¬ 
sults are based on the N2 method for the metallicity estimation. 
Although Andrews & Martini (2013) and Zahid et al. (2014b) 
adopted a projection parameter of a = 0.30, which minimizes 
the scatter in metallicities of their mass-metallicity relations, 
here we use a projection parameter of a = 0.32 for compari¬ 
son. The results change very little, even if we use a = 0.30. 


In the right panel of figure 6, we also show results from ear¬ 
lier works at ^ ^ 1.4 — 1.6 by Yabe et al. (2014) and Zahid 
et al. (2014b). All stellar masses and SFRs are based on the 
Salpeter IMF, and the metallicity calibration of the N2 method 
is assumed. 

The right panel of figure 6 shows that the metallicity of our 
sample is close to that of the local FMR by Mannucci et al. 
(2010), where the SFR is not corrected for the aperture effect, 
at higher /i. On the other hand, at lower p, the metallicity is 
larger than the local value by up to ^ 0.2 dex. Our result at 
2 ; ~ 1.4 is located between the FMRs with aperture corrected 
SFR at 2 ; 0.1 by Andrews & Martini (2013) and by Zahid et 

al. (2014b). 

4.3. Electron Density 

As presented in figure 5, we have significant detection of 
[S h]AA 67 16,6731 doublet in the stacked spectra. The obtained 
[S h]A6716/[S h]A 6731 line ratio ranges from ^ 1.0 to 1.8 
with a median value of 1.33. The corresponding electron den¬ 
sity ranges rie ^ 10 to 500 cm“^ with a median value of 80 
cm“^ assuming the electron temperature of 10,000 K (Shaw & 
Dufour 1995). In this range of the line ratio, the electron den¬ 
sity does not change significantly (up to a factor of 2) even 
if we assume the electron temperature of 100,000 K. The re¬ 
sulting line ratios, which are also summarized in table 2, are 
plotted against the stellar mass in figure 7. No clear trend be¬ 
tween the electron density and other physical parameters such 
as stellar mass and SFR can be seen. The range of the line ratio 
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Table 2. Emission line ratios, metallicity, N/0 in each stellar mass and SFR bin of our sample. 
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No. ] 


Bin 

N2 index® 

S2 index'' 

S2 ratio® 

12H-log(0/H)® 

log(N/0)® 

12+log(0/H)eo./ 

1 

-1.03±0.06 

-0.77±0.08 

0.88±0.36 

8.31±0.01 

-1.20±0.12 

8.13±0.01 

2 

-0.91±0.04 

-0.81 ±0.05 

1.98±0.55 

8.38±0.01 

-0.98±0.08 

8.23±0.01 

3 

-0.73±0.04 

-0.60±0.06 

1.71±0.53 

8.48±0.02 

-1.02±0.08 

8.37±0.02 

4 

-0.66±0.07 

-0.76±0.07 

1.38±0.46 

8.52±0.03 

-0.74±0.12 

8.43±0.03 

5 

-0.82±0.08 

-0.77±0.09 

1.46±0.71 

8.43±0.02 

-0.93±0.14 

8.30±0.02 

6 

-0.99±0.06 

-0.84±0.06 

0.93±0.29 

8.34±0.01 

-1.04±0.10 

8.15±0.01 

7 

-0.87±0.04 

-0.65±0.04 

1.51±0.30 

8.40±0.01 

-1.14±0.07 

8.24±0.01 

8 

-0.82±0.05 

-0.62±0.05 

1.35±0.29 

8.43±0.02 

-1.11±0.08 

8.28±0.02 

9 

-0.70±0.05 

-0.70±0.05 

1.35±0.34 

8.50±0.02 

-0.86±0.09 

8.38±0.02 

10 

-0.74±0.06 

-0.73±0.09 

1.18±0.52 

8.48±0.02 

-0.87±0.13 

8.34±0.02 

11 

-0.92±0.04 

-0.71 ±0.05 

1.71±0.39 

8.38±0.01 

-1.12±0.08 

8.19±0.01 

12 

-0.89±0.04 

-0.74±0.05 

1.75±0.47 

8.39±0.01 

-1.06±0.08 

8.21±0.01 

13 

-0.73±0.04 

-0.69±0.06 

1.07±0.30 

8.48±0.02 

-0.92±0.09 

8.34±0.02 

14 

-0.69±0.04 

-0.69±0.06 

1.78±0.49 

8.51 ±0.02 

-0.85±0.08 

8.37±0.02 

15 

-0.66±0.05 

-0.80±0.07 

1.32±0.45 

8.52±0.02 

-0.68±0.10 

8.39±0.02 

16 

-0.99±0.06 

-0.89±0.06 

1.12±0.33 

8.33±0.01 

-0.99±0.10 

8.12±0.01 

17 

-0.85±0.05 

-0.74±0.05 

1.22±0.30 

8.41±0.01 

-1.00±0.08 

8.22±0.01 

18 

-0.82±0.04 

-0.69±0.04 

1.29±0.26 

8.43±0.01 

-1.02±0.06 

8.25±0.01 

19 

-0.70±0.03 

-0.70±0.05 

1.45±0.33 

8.50±0.01 

-0.86±0.07 

8.35±0.01 

20 

-0.70±0.04 

-0.81 ±0.08 

1.17 ±0.44 

8.50±0.01 

-0.71±0.10 

8.35±0.01 

21 

-1.03±0.07 

-0.85±0.07 

1.40±0.45 

8.31±0.01 

-1.09±0.11 

8.06±0.01 

22 

-0.86±0.05 

-0.77±0.05 

1.62±0.37 

8.41±0.01 

-0.97±0.08 

8.19±0.01 

23 

-0.80±0.04 

-0.80±0.05 

1.42±0.32 

8.45±0.01 

-0.85±0.07 

8.25±0.01 

24 

-0.66±0.04 

-0.70±0.04 

1.16±0.21 

8.52±0.02 

-0.82±0.06 

8.35±0.02 

25 

-0.69±0.04 

-0.68±0.05 

1.21±0.26 

8.51±0.01 

-0.87±0.07 

8.33±0.01 


Note - The bin number corresponds to that in Table 1. Errors are estimated by using the bootstrap method, (a) log([N 
ll]A6583/Ha) (b) log([S ll]AA6716,6731/Ha) (c) [S ll]A6716/[S ll]A6731 (d) Oxygen abundance determined from the N2 index 
(e) N/O abundance ratio determined from the N2S2 index (=log([N ll]A6583/[S ll]AA6716,6731)). (f) Oxygen abundance after 

correcting the effect of N/O abundance ratio. 


of our sample is in good agreement with that of local galaxies 
selected from the SDSS DR8 catalogue (Aihara et al. 2011), 
suggesting that galaxies at z ^ 1.4 have similar electron den¬ 
sity to local galaxies. The resulting electron density at z ^ 1.4 
roughly agrees with the recent result (ue = 100 — 400 cm“^) 
of low mass galaxies at z ^ 2 by Masters et al. (2014). 

5. Discussion 

5.7. Nitrogen-to-Oxygen Abundance Ratio 

There is a concern that the metallicity estimation has sys¬ 
tematic errors (e.g., Nagao, Maiolino, & Marconi 2006; 
Kewley & Ellison 2008). The estimation based on the 
[N ll]A6583 /Ha line ratio depends on the ionization parame¬ 
ters (Denicolo, Terlevich, & Terlevich 2002; Kewley & Dopita 
2002). Recently, Masters et al. (2014) examine emission line 
galaxies at z 2, and find an offset from local galaxies on a 
[O lll]AA4959,5007/H/3 vs. [N ll] A6583 /Ha diagram, also re¬ 
ported in earlier studies (e.g., Kewley et al. 2013), but no clear 
offset on a [O iii]AA4959,5007/H/3 vs. [S ii]AA6716,6731/Ha 
diagram. They conclude that the significant offset on the 
[O lll]AA4959,5007/H/3 vs. [N ll]A6583/Ha diagram is caused 
by a nitrogen enhancement due to the high nitrogen-to-oxygen 
abundance ratio (hereafter N/O). Steidel et al. (2014) suggest 


that the higher N/O ratio of high-redshift galaxies can partly ex¬ 
plain a similar offset of high-redshift galaxies from local galax¬ 
ies on the BPT diagram. Hayashi et al. (2015) also show that 
the N/O of galaxies at z ^ 1.5 derived by using [N ll] A6583 and 
[O ll]A3727 lines is > 0.1, which is higher than local galaxies 
at a fixed metallicity. In this subsection, we will examine the 
N/O of galaxies at z 1.4 using emission lines measured from 
the composite spectra of our large NIR spectroscopic sample. 

Eigure 8 shows the distribution of [N ll]A6583/Ha vs. 
[S ll]AA6716,6731/Ha in each stellar mass and SER bin. Both 
line ratios are often used in diagnostics to distinguish between 
star-forming galaxies (SEGs) and AGNs (e.g., Osterbrock 
1989). Empirical separation lines between SEGs and AGNs 
on [N II]A6583/Ha vs. [S ll]AA6716,6731/Ha diagram, which 
are indicated in figure 8, have been proposed in previous stud¬ 
ies (Lamareille et al. 2009, and Lara-Lopez et al. 2010a). Our 
data points are located in the region of the SEGs, which is suffi¬ 
ciently far from the AGN region, implying that the contamina¬ 
tion of AGN to our sample would be very low. Note that broad 
line AGNs with the line width of > 1000 km s“^ are excluded 
from our sample beforehand. 

Eigure 8 shows that the distribution of our sample dif¬ 
fers from that of SDSS normal SEGs; In our sample, both 
[N II]A6583/Ha and [S ll]AA6716,6731/Ha are lower than 
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those of main part of the SDSS sample. In figure 8, our sam¬ 
ple is compared to emission line galaxies at z ^ 0.1 to ^ 0.3 
called Green Pea galaxies (GPs; Cardamone et al. 2009). Here 
we use the GP sample by Cardamone et al. (2009) with the 
emission line fluxes summarized by Hawley (2012). The GPs 
are generally located in the region with lower [N II]A6583 /Hq: 
and [S ll]AA6716,6731/Ha than the SDSS sample. In the high 
[N ii]A6583/Ha and [S ii]AA6716,6731/Ha end, they overlap 
with our sample. 


in 



Fig. 7. The distribution of [S II]A6716/[S II]A6731 line ratio against 
the stellar mass. Symbols are color-coded in the same manner as fig¬ 
ure 6. The contour shows the distribution of the local SDSS galax¬ 
ies. The horizontal dotted lines are the corresponding electron den¬ 
sities. Here, we assume the electron temperature of 10,000 K. Gray 
shaded region indicates the results at 2 : 2 by Masters et al. (2014). 


The correlation between [N ll]A6583 to [S ll]AA6716,6731 
line ratio and N/O has been investigated in local galaxies and 
empirical calibrations have been proposed (Perez-Montero & 
Contini 2009, Amorm, Perez-Montero & Vilchez 2010, and 
Hawley 2012). In this work, we estimate the N/O ratio for our 
sample at z ~ 1.4 by using the N2S2 index calibrated by Perez- 
Montero & Contini (2009), where the N2S2 index is defined as 
follows: 


N2S2 = log 


/ /([Nii]A6583) A 

V([Sii]AA6716,6731 )) ' 


( 2 ) 


Perez-Montero & Contini (2009) investigate the relation be¬ 
tween the N2S2 index and the N/O of local galaxies. Although 
the scatter of the relation is moderately large (cr ^ 0.3 dex), 
they present the least-squares bisector linear fit: 


log(N/0) = 1.26 X N2S2 - 0.86. 


(3) 


In figure 9, the N/O measured in each stellar mass and SFR 
bin of our sample galaxies at z ^ 1.4 is presented as a function 
of metallicity and stellar mass, which is summarized in table 2. 
The mean (median) log(N/0) is —0.95 (—0.98). The N/O ratio 
increases with increasing both metallicity and stellar mass. The 



Fig. 8. [N II]A6583/Hq: versus [S II]AA6716,6731 /Hq!. Our samples 
with different SFR groups are indicated as filled circles with the same 
color-coding as figure 6. The larger symbols show the higher stel¬ 
lar masses. The distribution of the SFGs and AGNs taken from 
the SDSS DR8 catalogue is shown by red and blue contour, respec¬ 
tively. The empirical lines separating SFGs and AGNs by Lamareille 
et al. (2009) and Lara-Lopez et al. (2010a) are shown by dotted 
and dashed line, respectively. The Green Pea galaxies (Cardamone 
et al. 2009, and Hawley 2012) are indicated by open circles. 


N/O ratio also appears to increase with increasing SFR at a 
fixed metallicity and stellar mass, although the trend is only 
marginal. 

The obtained N/O of our sample at z ^ 1.4 is compared 
with those of local galaxies in previous observations (Perez- 
Montero & Contini 2009, Pilyugin et al. 2012, Berg et al. 2012, 
Andrews & Martini 2013, and Dopita et al. 2013), where the 
stellar masses are scaled to those derived assuming the Salpeter 
IMF. The N/O at z ~ 1.4 is generally higher than the local val¬ 
ues at a fixed metallicity but almost comparable to those at a 
fixed stellar mass, which may be probably due to the difference 
of the mass-metallicity relation between local and high-redshift 
galaxies. The resulting N/O of galaxies at z ~ 1.4 well agrees 
with those of massive and metal-rich part of the GPs sample 
by Hawley (2012). It is worth noting that the deviation of N/O 
from local galaxies is significantly larger than the intrinsic un¬ 
certainty of the N/O calibration using the N2S2 index. 

Queyrel et al. (2009) estimated the N/O of four star-forming 
galaxies at z = 1.3 - 1.5 by using the same N2S2 index. The 
obtained log(N/0) ranges from —0.5 to —1.1 in the metallicity 
range of 12H-log(0/H) = 8.4 - 8.6 and the stellar mass range 
of 10^° - 10^^ Mq, which are almost similar to those of our 
sample. More recently, the similarly high N/O of log(N/0) ~ 
— 1.0 is also reported by using galaxy samples at z ~ 2 (Masters 
et al. 2014, Steidel et al. 2014, and Hayashi et al. 2015). 

One of the natural explanations of the variation in the N/O 
may be a time delay between the release of oxygen and that 
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Fig. 9. (Top) N/0 vs. 12+log(0/H). The symbols are the same as 
in figure 6. Galaxies with the larger stellar masses are indicated by 
the larger symbols. The thick solid, dashed, dot-dashed, and dot¬ 
ted lines are results for local galaxies by Andrews & Martini (2013), 
Berg et al. (2012), Dopita et al. (2013), and Pilyugin et al. (2012), 
respectively. The open squares are results at z = 1.3 — 1.5 by 
Queyrel et al. (2009). (Bottom) N/0 versus stellar mass. Thick 
solid and dashed line are results by Perez-Montero & Contini (2009) 
and Andrews & Martini (2013), respectively. The vertical and 
horizontal dashed line indicate the solar abundance and the so- 
lai' abundance ratio, respectively (Asplund et al. 2009). The open 
squares are results at z = 1.3 — 1.5 by Queyrel et al. (2009). 


of nitrogen into the ISM; the two elements are formed in stars 
of different masses with different lifetimes (e.g., Edmunds & 
Pagel 1978, Garnett 1990, Kobulnicky & Skillman 1998, and 
van Zee, Salzer, & Haynes 1998). In this scenario, after a sin¬ 
gle starburst, the N/O decreases and O/H increases due to the 
release of oxygen formed in massive, short-lived stars with the 
time-scale of a few tens of Myr. Then, the delayed release of 
nitrogen produced in intermediate mass, longer-lived stars in¬ 
creases N/O over longer time period of a few hundred Myr, 
while the O/H does not change significantly. If there are mul¬ 
tiple starburst events, this cycle is repeated and both N/O and 
O/H increase gradually (e.g., Garnett 1990, and Contini et al. 
2002). The variance in the distribution of the N/O of local 
galaxies could be explained by the combination of this cycle. 
The high-redshift galaxies with high N/O may correspond to 
the upper locus of the distribution of local galaxies, where they 
are in the late quiescent phase of the nitrogen release. It may 
be unlikely, however, that all of the high-redshift galaxies are 
located in the quiescent phase. 

IMF can change the chemical abundance ratio (Kobulnicky 
& Skillman 1998). Since the nitrogen predominantly originates 
from low- or intermediate-mass stars, the N/O increases if the 
IMF has a steep slope or bottom-heavy shape. In recent stud¬ 
ies, the evolution of the IMF with redshift has been suggested 
(e.g., van Dokkum 2008). The slope of the massive end of the 
IMF at high redshift may depends on the physical properties 
of galaxies such as luminosity (Hoversten & Glazebrook 2008) 


and SFR (Gunawardhana et al. 2011). The trend of the higher 
N/O at high redshift or high SFR regime, however, does not 
agree with the recent suggestions of flatter or top-heavy IMF. 

The N/O enhancement by the nitrogen-rich (WN) Wolf- 
Rayet (WR) stars (e.g., Crowther 2007) could be another cause 
of the higher N/O at high redshift (e.g., Kobulnicky & Skillman 
1998, Henry, Edmunds, & Koppen 2000, and Brinchmann, 
Kunth, & Durret 2008). Increasing contribution of WR stars 
is expected in flatter IMF (e.g., Schaerer, Contini, & Kunth 
1999). In some studies of nitrogen enhanced galaxies, the 
spectroscopic features originated from WR stars have been re¬ 
ported (Pilyugin et al. 2012, Brinchmann, Kunth, & Durret 
2008, and Amorln et al. 2012). Perez-Montero et al. (2011), 
however, suggest no positive evidence of local pollution of ni¬ 
trogen enriched gas from the WR stars in local blue compact 
dwarf galaxies with the integral held spectroscopy. 

Another possible explanation is the dilution effect by the in¬ 
fall of metal poor gas (e.g., Edmunds 1990). In this scenario, 
galaxies were originally more metal-rich and the distribution 
was close to the metal-rich part of local galaxies, where N/O in¬ 
creases with increasing metallicity because the secondary pro¬ 
duction of nitrogen is dominated (Koppen & Hensler 2005, van 
Zee & Haynes 2006, and Masters et al. 2014). Then, the in¬ 
fall of the metal poor gas causes the decrease of the metal¬ 
licity while N/O ratio does not change. In recent studies, es¬ 
pecially in the theoretical context, gas inflow plays an impor¬ 
tant role in the galaxy growth. According to numerical models 
by Henry, Edmunds, & Koppen (2000), however, the time re¬ 
quired to approach the line of secondary production regime on 
12H-log(0/H) vs. N/O diagram is, at least, > 2 Gyr, and thus, it 
may be unlikely that all of our sample galaxies at ^ 1.4 were 
initially located in the secondary production regime of nitrogen 
(see discussion by Masters et al. 2014). 

The differential outflow of oxygen-rich gas could be the ori¬ 
gin of the higher N/O at high redshift. Recently, various studies 
suggest that the gas outflow is ubiquitous in star-forming galax¬ 
ies at high redshift (e.g., Weiner et al. 2009, Steidel et al. 2010, 
and Yabe et al. 2015). Outflowing gas originated from super¬ 
nova (SN) wind could be rich in oxygen relative to nitrogen be¬ 
cause oxygen is generally formed in high-mass and short-lived 
stars while nitrogen is supplied from intermediate-mass and 
longer-lived stars. Although the gas outflow with well-mixed 
material does not change the abundance pattern of galaxies, 
the selective outflow of oxygen-rich gas driven by SN could 
cause the nitrogen enhancement (van Zee & Haynes 2006). 
The direct observational evidence of the oxygen-rich gas out¬ 
flow, however, would be desirable for further arguments. 

5.2. The Mass-Metallicity Relation and the Fundamental 
Metallicity Relation Revisited 

In the previous subsection, we show the higher N/O in high- 
redshift galaxies than local galaxies at a fixed oxygen abun¬ 
dance. This nitrogen enhancement causes a systematic effect 
on the metallicity estimation with the N2 method based on the 
[N ll]A6583/Ha flux ratio; if the N/O is higher than the local 
value, the metallicity derived from the N2 method is overesti¬ 
mated. Perez-Montero & Contini (2009) obtain the alternative 
N2 method taking into account the N/O effect: 
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Fig. 10. The mass-metallicity relation (left) and the fundamental metallicity relation (right) based on the recalculated metallicity of our sample 
and those by Yabe et al. (2014) and Zahid et al. (2014b) corrected for the N/0 enhancement. The symbols are all the same as figure 6. 


12+log(0/H) = 0.79 X N2 - 0.56 x log(N/0) + 8.41, (4) 

where N2 = log([N ll]A6583/Ha). This empirical calibration 
shows that the metallicity decreases with increasing N/O at a 
fixed N2 value. We recalculate the metallicity of our sample at 
z ^ 1.4 by using the above calibration taking into consideration 
the N/O enhancement. We use the mean value of log(N/0) 
= —0.95 described in subsection 5.1 for all of our sample. 

The metallicity corrected for the differing N/O is generally 
lower than the original metallicity by 0.1 - 0.2 dex, which is 
also summarized in table 2. The left panel of figure 10 shows 
the distribution of the recalculated metallicity of our sample 
and previous results at z ^ 1.4 — 1.6 by Yabe et al. (2014) and 
Zahid et al. (2014b) against the stellar mass. Here, the N/O ra¬ 
tio of previous studies is assumed to be the same value as our 
sample, and the metallicities are corrected in the same way as 
described above. The recalculated metallicities of our result are 
generally lower than that of the result at z ~ 0.1 at a fixed stel¬ 
lar mass by up to ^ 0.4 dex. The right panel of figure 10 shows 
the distribution of the corrected metallicity of our sample on 
the FMR diagram. Although there exists a discrepancy of up to 
~ 0.3 dex in the larger p, the result with the corrected metallic¬ 
ity agrees with the extrapolation of the local FMR by Mannucci 
et al. (2010), where the SFR is not corrected for the aperture ef¬ 
fect, in the other part. On the other hand, the resulting metallic¬ 
ity corrected for the N/O enhancement is lower than the results 
at z ~ 0.1 with aperture corrected SFR by Andrews & Martini 
(2013) and Zahid et al. (2014b) at fixed fj,, implying that the 
FMR may slightly evolve from z ^ 1.4toz^0.1. 

6. Summary 

We present the results from a large near-infrared spectro¬ 
scopic survey with Subaru/FMOS {FastSound) consisting of ^ 
4,000 galaxies at z ~ 1.4 with significant Ha detection. We 
measure the gas-phase metallicity from the [N ll]A6583/Ha 


emission line ratio of the composite spectra in various stel¬ 
lar mass and star-formation rate bins. The resulting mass- 
metallicity relation generally agrees with previous studies ob¬ 
tained in a similar redshift range to that of our sample. The 
metallicity of our sample derived with the N2 method is lower 
than that at z ^ 0.1 by up to 0.2 dex at a fixed stellar mass. 
No clear dependence of the mass-metallicity relation with star- 
formation rate is found. The metallicity of our sample is sig¬ 
nificantly larger than the FMR at z ^ 0.1 derived by Mannucci 
et al. (2010) especially in lower p. The resulting metallicity 
roughly agrees with the FMRs at z ^ 0.1 by Andrews & Martini 
(2013) and Zahid et al. (2014b), where the SFR is corrected 
for the fiber aperture effect and the metallicity is estimated 
with the N2 method. We detect significant [S ll]AA6716,6731 
emission lines from the composite spectra. The electron den¬ 
sity estimated from the [S ll]AA6716,6731 line ratio ranges 
from 10 - 500 cm“^, which generally agrees with that of local 
galaxies. On the other hand, the distribution of our sample on 
[N II]A6583/Ha vs. [S ll]AA6716,6731/Ha is different from 
that found locally, but is similar to that of lower redshift emis¬ 
sion line galaxies such as Green Pea galaxies. We estimate the 
nitrogen-to-oxygen abundance ratio (N/O) from the N2S2 in¬ 
dex, and find that the N/O in galaxies at z ~ 1.4 is significantly 
higher than the local values at a fixed metallicity and stellar 
mass. The metallicity of our sample at z ^ 1.4 recalculated 
with this N/O enhancement taken into account is lower than 
the original value by 0.1 - 0.2 dex and is lower than the local 
value by up to 0.4 dex at a fixed stellar mass. The recalculated 
metallicity roughly agrees with the local FMR by Mannucci et 
al. (2010) but is lower than the FMRs with the aperture cor¬ 
rected SFRs by Andrews & Martini (2013) and Zahid et al. 
(2014b). 
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